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A  nevr  technique  for  full-order  and  reduced-order  observer 
design  which  is  based  on  the  minimization  of  an  integral 
square-error  performance  measure  is  developed  for  the  inte- 
gral of  time-multiplied-square-error  performance  measure. 
The  technique  is  extended  to  include  the  design  of  full- 
order  and  reduced-order  combined  observer/controllers.   The 
technique  is  well  suited  to  the  design  of  higher  order, 
multiple -input ,  multiple-output  systems  where  traditional 
methods  become  difficult  to  apply. 
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I.   INTRODUCTION 

In  feedback  control  systems,  the  plant  control  is  often 
a  function  of  the  current  plant  state  vector.   A  knowledge 
of  the  state  vector  is  therefore  required  for  the  controller. 
In  most  systems,  however,  the  entire  state  vector  is  not 
available  through  direct  measurement.   It  is  therefore 
necessary  to  estimate  the  unknown  states.   An  observer  is  a 
device  which  can  be  used  to  reconstruct  all  or  part  of  the 
state  vector  [1] . 

Traditionally,  observers  and  controllers  have  been 
designed  separately.   Time  and  frequency  domain  specifica- 
tions are  satisfied  through  the  selection  of  controller 
feedback  gains  as  if  all  of  the  plant  states  are  available. 
The  observer  is  designed  to  estimate  the  plant  states  by 
proper  selection  of  the  observer  gain  matrix.   The  observer 
output  is  used  to  supply  estimates  of  the  missing  states  to 
the  controller  [2]. 

In  this  thesis,  after  presenting  the  necessary  background 
material,  a  method  for  designing  observers  which  minimizes 
a  quadratic  performance  measure  is  presented.   The  method  is 
also  extended  to  the  simultaneous  design  of  the  observer  and 
the  controller.   For  simplicity,  only  linear,  time- invariant 
systems  are  considered.   The  design  method  relies  on  the  use 
of  a  digital  computer  algorithm,  however,  the  computations 
are  performed  off-line.   Numerical  examples  are  presented 
in  order  to  illustrate  the  new  technique. 


II.   MATHEMATICAL  BACKGROUND 

A.   TRADITIONAL  DESIGN  TECHNIQUES 

1"  Vi 

The  dynamic  characteristics  o£  an  n   -order,  linear, 
time- invariant ,  free  system  can  be  described  by  the  linear, 
first  order  differential  equations 

s(t)  =  HsCt)  (1) 

where  s_(t)  is  an  nxl  vector  and  H  is  a  constant  nxn  matrix. 
The  response  of  the  system  to  any  non-zero  initial  condition 
is  determined  by  the  elements  of  the  matrix  H  which  are,  in 
turn,  a  function  of  the  system  configuration  and  physical 
constants.   By  adjusting  system  parameters,  a  suitable  H  ma- 
trix is  realized  which  causes  the  system  to  meet  certain 
system  performance  specifications. 

Traditionally,  performance  is  specified  as  a  consistent 
set  of  time  and  frequency  domain  criteria.   Typical  speci- 
fications are  given  in  terms  of  rise  time,  settling  time, 
percent  overshoot,  steady-state  error,  gain  margin,  phase 
margin,  and  bandwidth.   These  criteria  are  then  translated 
into  a  control  strategy  through  the  use  of  some  relatively 
simple  mathematical  tools;  generally  by  trial-and-error 
techniques.   Appropriate  assumptions  are  usually  made  to 
attempt  to  approximate  complex,  higher-order  systems  by 
simpler,  lower-order  models  that  can  be  handled  with  greater 
ease . 

This  approach,  although  empirical  in  nature,  and  rely- 
ing heavily  on  the  ingenuity  and  system  know-how  of  the 
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designer,  is  relatively  straightforward  for  simple,  single- 
input,  single-output  systems.   Many  modern  systems  are 
highly  complex,  however,  with  multiple  inputs  and  multiple 
outputs.   These  systems  do  not  readily  lend  themselves  to 
the  mathematical  techniques  alluded  to  earlier  and  often 
performance  is  not  easily  specified  in  terms  of  the  tradi- 
tional criteria.   It  is  therefore  desirable  to  have  an  al- 
ternative method  for  prescribing  desired  system  performance. 

B.   THE  QUADRATIC  PERFORMANCE  MEASURE 

One  such  index  of  system  performance  is  the  quadratic 
performance  measure 

00 

J   =   J    t^^(t)qs(t)dt  (2) 

o 

which  assigns  a  positive  scalar  value  to  the  infinite  time 
integral  of  the  time-multiplied  quadratic  form  in  £(t)  and 
Q^.   The  elements  of  the  nxn  matrix  Q^   are  specified  to  deter- 
mine the  relative  contributions  to  J  of  the  elements  of  the 
vector  s_(t)  .   The  value  of  the  time-multiplier  exponent,  r, 
establishes  the  relative  contributions  to  J  of  the  earlier 
and  later  values  of  the  elements  of  s_Ct)  .   Interpreting  J 
as  a  cost  associated  with  a  particular  system  configuration, 
it  becomes  the  task  of  the  designer  to  find  values  for  the 
elements  of  H  which  are  at  his  disposal  that  yield  the  lowest 
value  for  J.   This  choice  for  H  is  said  to  be  optimal  since 
it  minimizes  the  cost. 

Examples  of  systems  to  which  a  performance  measure  of 
the  form  of  (2)  can  be  applied  are  the  observer  error  state 


equations 


e(t)  =  FeCt) 


C3) 


and  the  state-error  differential  equations  of  a  combined 
plant -observer /controller 


xCt) 
e(t) 


A+BK  I  BK 


x(t) 
e(t) 


(4) 


Equations  (3)  and  (4)  will  be  developed  in  detail  in  Sections 
III  and  IV,  respectively. 

C.   EVALUATING  THE  QUADRATIC  PERFORMANCE  MEASURE 

In  Ref.  3  it  is  shown  that  if  (1)  is  a  stable,   linear 
system  and  if  ^   is  any  positive  definite  matrix,  the  per- 
formance measure  (2)  can  be  expressed  as 


J  =  (-D^^^r:  J(0)V^^^s(0) 


C5) 


The  real  symmetric  matrix  V  ^^  is  the  unique  solution  of 
the  simultaneous  linear,  algebraic  matrix  equations 


H^V  ^1  +  V  ^^H  =  V 
s  +  1   — s  +  1—   — s 


s  =  1,2, 


(6) 


It  can  be  shown  [4]  that  if  H  is  a  stability  matrix  and  if 
^   is  positive  definite  then  the  matrix  V^  is  also  positive 

definite.   The  argument  can  easily  be  extended  to  show  that 

r+ 1 
(-1)    rl V  ^,  is  also  positive  definite. 


The  system  (1)  is  stable  if  H  is  a  stability  matrix; 
that  is  if,  and  only  if,  all  eigenvalues  of  H  have  negative 
real  parts . 
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The  details  of  the  derivation  of  (5)  and  (6)  are  given 
in  Appendix  A.   The  solution  to  (6)  can  be  obtained  by 

iteratively  solving  the  r+1  sets  of  n(n+l)/2  simultaneous 

2 
linear  equations   or  by  a  technique  suggested  in  Ref.  3. 

The  details  of  this  technique  are  given  in  Appendix  B. 

Evaluation  of  the  performance  measure  (2)  using  equa- 
tions (5)  and  (6)  eliminates  the  need  for  determining  the 
explicit  time  response  of  the  system  but  requires  knowledge 
of  the  initial  condition,  s_(0)  .   Since  a  priori   knowledge 
of  this  vector  is  not  generally  available,  the  question 
arises  as  to  how  to  pick  s_CO)  in  order  to  arrive  at  a  mean- 
ingful value  for  J. 

One  technique,  suggested  in  Ref.  5,  is  to  determine  the 
average  value  of  J  over  a  linearly  independent  set  of  ini- 
tial condition  vectors.   This  is  equivalent  to  assuming 
that  the  initial  condition  vector  is  uniformly  distributed 
on  the  surface  of  a  unit  hypersphere  and  results  in  a  de- 
sign which  is  optimal  in  an  average  sense.   Another  tech- 
nique, suggested  in  Ref.  6,  is  to  determine  the  largest 
possible  J  for  every  s_(0)  within  or  on  a  hypersphere  of 
radius  p  and  select  the  configuration  which  yields  the 
smallest  of  these  maximum  values.   This  results  in  a  design 
which  is  optimal  in  a  worst  case  sense.   For  the  latter 
scheme,  it  is  shown  in  Appendix  C  that  for  a  given  design 
parameter  set  {P},  the  largest  value  for  J  is  given  by  J. 
where 


2 
Since  the  matrices  V,,  V7.---,V,+1  are  each  symmetric, 

only  n(n+l)/2  elements  of  each  need  to  be  solved  for. 
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and  where  ^^^^^.f*]  represents  the  largest  eigenvalue  of  [•]. 
Here,  the  functional  dependence  of  J,  upon  the  design 
parameter  set  {P}  is  stressed. 

D.   MINIMIZING  THE  QUADRATIC  PERFORMANCE  MEASURE 

To  this  point  we  have  assumed  that  the  design  parameter 
set  {P}  is  unbounded.   In  general,  however,  we  must  allow 
for  the  possibility  of  a  constrained  {P}.   There  may  be 
explicit  as  well  as  implicit  constraints.   The  explicit  con- 
straints are  those  which  are  the  result  of  the  physical 
limitations  on  the  system  components  while  the  implicit 
constraints  are  those  which  result  from  restrictions  placed 
on  system  performance  by  the  designer 

In  either  case,  the  constraining  relationships  must  be 
included  in  the  problem  formulation.   One  means  of  accom- 
plishing this  is  to  add  a  large  penalty  to  the  cost  of  a 
given  configuration  if  that  configuration  requires  the  vio- 
lation of  any  of  the  constraints.   With  this  in  mind  we  de- 
fine the  optimal  parameter  set  {P  }  as  that  parameter  set 

which  minimizes  J^({P})  and  the  minimum  value  of  J,({P*}) 

* 
as  J  ;  that  is 

J*=  J,({P*})  =  min  Jt({P}). 

{P}   ^ 

■z 

Therefore 


3c-      ,  - 

bmce  r.  is  a  constant  in  any  given  design  specification, 

It  docs  not  contribute  to  the  optimization  of  the  parameter 
set  and  can  therefore  be  deleted  from  the  performance  measure. 
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This  indicates  that  the  following  iterative  procedure 
can  be  used  to  find  {P  }  once  r  and  Q^   have  been  specified. 

1.  Select  an  initial  estimate  for  {P}  insuring  that 
any  constraints  imposed  are  satisfied  and  that  H  is  a  sta- 
bility matrix. 

2.  Using  equation  (6)  or  the  method  described  in  Ap- 
pendix B,  find 

3.  Evaluate 

4.  Through  some  suitable  minimizing  search  routine 
select  a  new  set  of  parameters  {P}.   If  the  constraints  are 
not  satisfied,  add  a  penalty  to  the  cost  and  select  a  new 
parameter  set. 

5.  Repeat  steps  two  through  four  until  some  criterion 
for  ending  the  search  is  satisfied. 


4 
For  the  examples  used  in  this  thesis,  the  subroutine 

"DIRECT"  (NPS  Computer  Facility)  tliat  performs  the  pattern 

search  method  of  Ilooke  and  Jeeves  was  used. 
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III.   OBSERVER  DESIGN 

In  the  design  of  feedback  control  systems  it  is  often 
desirable  to  have  information  about  plant  states  which  is 
not  normally  available  through  direct  measurement.   It  is 
sometimes  possible  to  construct  a  system  called  an  observer 
which  can  estimate  a  linear  transformation  of  the  plant 
state  and,  if  the  transformation  is  invertible,  it  is  then 
possible  to  reconstruct  the  plant  state  vector  itself. 

Traditionally,  observer  design  has  been  accomplished 
through  the  use  of  techniques  such  as  those  described  in 
Section  IIA.   In  this  section,  techniques  for  the  design  of 
full-order  and  reduced-order  observers  by  using  the  quad- 
ratic performance  measure  as  the  design  specification  are 
presented. 

A.   FULL-ORDER  OBSERVER  DESIGN 

We  consider  the  time-invariant,  linear  system 

x(t)  =  Ax(t)  +  Bu(t) 

(9) 
X(t)  =  Cx(t) 

where  x(t)  is  the  nxl  state  vector,  )^(t)  is  the  mxl  output 
vector,  u(t)  is  an  Ixl    control  vector,  A  is  an  nxn  constant 
matrix,  B  is  an  nx£  constant  matrix,  and  C  is  an  mxn  con- 
stant matrix.   We  consider  also  the  full-order  identity 
observer  shown  in  Fig.  1.   The  observer  state  equations  as 
developed  in  Appendix  D  are  given  by 

z(t)  =  (A-GC)£(t)  +  G^(t)  +  Bu(t)        (10) 
14 


xCt) 


x(t) 


I(t) 


u(t) 


Figure  1.   Full-Order  Observer 
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where  z^(t)  is  the  nxl  observer  state  vector,  G  is  the  nxm 
observer  gain  matrix  (in  this  case  the  elements  o£  G_   con- 
stitute the  design  parameter  set  {P})  ,  and  where  it  is  as- 
sumed that  (C^jA)  is  observable. 
Defining  the  error  state  as 

A 
eCt)  =  z_(t)  -  x(t) 

and  differentiating  with  respect  to  time  yields 

e(t)  =  £Ct)  -  xCt)  . 

Substituting  (9)  and  (10)  for  x(t)  and  z_(t)  ,  respectively, 
and  solving  for  ^(t)  yields 

e(t)  =  (A-GC)e(t). 

We  define  the  matrix  F  as 

A 
F  =  (A-GC)  (11) 

which  yields 

i(t)  =  Fe(t).  (12) 

We  define  the  quadratic  performance  measure  as 

oo 

J   -    I      t^e'^(t)Se(t)dt.  (13) 

o 

The  system  (12)  and  the  performance  measure  (13)  have  the 
same  form  as  (1)  and  (2),  respectively.   Equation  (13)  can 
therefore  be  rewritten  as 

J  =  (-l)^^^r:e'^(0)V^^^e(0) 


where,  by  using  (6), 
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F  V  ^1  +  V  ^iF  =  V 
s  +  1   — s  +  1—   — s 


s  =  1,2, ... ,r 


Lli  *  liL      =  -2.- 

Furthermore,  the  algorithm  described  in  Section  IID  can  be 
applied  to  find  G^  ,  the  observer  gain  matrix  which  is  opti 
mal  vrith  respect  to  the  performance  measure  (13)  . 


B.   REDUCED-ORDER  OBSERVER  DESIGN 

We  again  consider  the  time- invariant ,  linear  system 
given  by 

x(t)  =  AxCt)  +  Bu(t) 

jrCt)  =  Cx(t). 


(14) 


We  define  the  vector  x(t)  as 

w(t) 


x(t)  ^ 


where  w(t)  is  an  (n-m)xl  vector  given  by 


w(t)  =  Dx(t) 


(15) 


and  where  D  is  an  (n-m)xn  matrix,  arbitrarily  selected  so 

that 

-1 
C 

D 

exists.   By  proper  selection  of  the  state  variables  C  will 
have  the  form 

C  =  [   I   :    0_  ]  . 

A  notationally  convenient  choice  for  D  is 

D  =  [  ^  :  ^  ]  . 
17 


We  can  find  w(t) ,  the  estimate  for  wCt)  by  constructing  the 
reduced- order  observer  as  shown  in  Fig.  2.   For  the  above 
choice  for  C  and  D,  the  observer  state  equations  developed 
in  Appendix  D  can  be  written  as 


+  CB2-GB^)u(t) 


where 


w(t)  =  ^(t)  +  Gx(t) 


(16) 


(17) 


and  where  the  matrices  A  and  B  are  partitioned  as 


A  = 


tII.    .t12 
^21  ^  -22 


B  = 


^1 

i2 


Here,  A^  ^  is  an  mxm  matrix,  A^  2  is  an  mx(n-m)  matrix,  A,-. 
is  an  (n-m)xm  matrix,  A-^  is  an  (n-m)x(n-m)  matrix,  B^  is 
an  mxil  matrix,  and  B^  is  an  (n-m)x£  matrix.   Therefore, 
jc_(t)  ,  the  estimate  for  ^(t)  (which  is,  for  the  case 
C  =  [  1  '  0^  ]  and  D  =  [  0^  :  1  ]  ,  the  estimate  for  x(t)), 
becomes 


x(t)  = 


zCt) 
z(t)  +  Gx(t)_ 

Defining  the  error  in  the  estimate  w(t)  as 

e(t)  ^  w(t)  -  w(t) 

and  substituting  (15)  and  (17)  for  w(t)  and  w(t) ,  respec 
tively,  yields 

e(t)  =  z(t)  +  (GC-D)x(t). 
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xCt) 


xCt) 


iCt) 


» 


A2rGA^lMA22-GA^2^G 


u(t) 


Figure  2,   Reduced-Order  Observer, 
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Differentiating  with  respect  to  time,  substituting  (14)  and 
(16)  for  xCt)  and  z^(t)  ,  respectively,  and  solving  for  £(t) 
yields 

e(t)  =  (A22-GA^2^^'^^^  *    ^^ZT^ll  *  ^t^zz'^lZ^^'^^^^^ 

+  (B2-GB^)u(t)  +  (GC-D) CAx(t)  +  BuCt)). 

We  note   that    for  C   =    [  I.  :    0.  ]    and  D   =    [    0_   :  X  ] 

GCAx(t)  =  [GA^j^    :    GA^2]2iCt) 

-DAx(t)  =  [-A2^    :    -A22]xCt) 

GCBu(t)  =  GB^uCt) 

-DBu(t)  =  -B2u(t). 

We  note     also   that    for  y_it')    =   Cx(t) 

A2ilCt)    =    [A2-^    :    0]x(t) 
-GA^llCt)    =    [-GA^^    :    0]x(t). 

These  relationships,  when  substituted  into  the  above  expres- 
sion  for  e  (t)  ,  yield 

e(t)  =  CA22-GA^2^^^^^  "  "^^2  2 '^^12^^'^''^^^ 

-[  0  :  (A22-GA^2^1^'^^^- 

We  note  that  the  last  term  in  this  expression  can  be  written 
as 

-  [  0  :  (A22-GA^2^1^'^^^  ^  '  (A22 -GA^2^  ^21^)  . 

Substituting  this  relationship  into  the  expression  for  £(t) 
and  combining  terms  yields 


eCt)  =  (A22-GA^2Hz(t)  +  (GC-D)x(t)) 
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where 

zCt)  +  CGC-D)xCt)  =  eCt). 

Therefore  e_it')    becomes 

l£  we  define 

F  ^  CA22-GA^2^  ^^^^ 

then 

e(t)  =  FeCt).  '       (19) 

Defining  the  quadratic  performance  measure  as 


CO 

J  =  I  t^e'^(t)gieCt)dt 


(20) 


we  see  that  the  system  (19)  and  the  performance  measure  (20) 
have  the  same  form  as  (1)  and  (2),  respectively,  and  that 
(20)  can  therefore  be  written  as 

J  =  (-l)^^^r:J(0)V^^^e(0). 
By  using  (6)  we  have 

fv^.l  :  V^.iF  =  V3        s  =  l,2,...,r 

Once  again,  the  algorithm  described  in  Section  IID  can  be 
used  to  find  G  ,  the  observer  gain  matrix  which  is  optimal 
with  respect  to  the  performance  measure  (20) . 

C.   EXAMPLE  OBSERVER  DESIGN  PROBLEM 

To  illustrate  the  application  of  the  algorithm  to  ob- 
server design  we  consider  the  plant  whose  transfer  function 
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is  given  by 


YfS 


(Sy  "  S3  +  12S2  +  47S  +  6 


0 


The  state  and  output  equations  for  this  plant  can  be  written 
as 


x(t)    = 


0  10 

0  0  1 

-60      -47      -12 


x(t)    + 


0 
0 
1 


u(t) 


)r(t)    =    [100    ]xCt). 


Here, 


A  = 


0  10 

0  0  1 

•60      -47      -12 


B   = 


0 
0 

1 


,    C    =    [    1    0    0    ] 


It  can  be  easily  shown  that  (C,A)  is  completely  observ- 
able.  A  second-order  observer  with  arbitrary  eigenvalues 
can  therefore  be  constructed  to  provide  estimates  for  the 
unmeasured  states  x- (t)  and  x_(t).   These  estimates,  x„(t) 
and  x_(t),  are  given  by 

x^Ct) 


X3(t) 


=  w(t)  =  ^(t)  +  GxCt) 


We  require  an  initial  estimate  for  the  observer  gain 
matrix  G_  which  yields  a  stable  observer.  The  observer  F 
matrix  is  given  by  (18)  as 


A 


22 


GA 


12 


(18) 


Here 


22 


-22 

0 

-47 

1 
-12 

Therefore 

,   A^2  =  [  1  0  ]  »    £  = 


gi 


1^2. 


F  = 


-gl    1 
■47-g2   -12 


where  the  observer  eigenvalues  are  the  eigenvalues  of  F  and 
are  given  by  the  solution  to 

X^  +  C12+g^)X  +  C12gi+g2+47)  =  0. 

We  arbitrarily  select  the  initial  estimate  for  the  ob- 
server eigenvalues  to  be  -7,  -2.  The  characteristic  equa- 
tion is  therefore 

(X+7) CX+2)  =  0 

or 

X^  +  9X  +  14  =  0. 

Equating  coefficients  of  the  like  powers  of  X  and  solving 
for  the  G  matrix  yields 


G  = 


•3 
3 


This  choice  for  G  must  also  satisfy  all  design  constraints 
A  parameter  in  observer  design  which  must  be  limited  is  the 
observer  speed  since,  as  the  observer  gains  are  adjusted  to 
cause  the  errors  to  die  out  more  rapidly,  the  observer  be- 
comes more  sensitive  to  noise.   One  measure  of  the  observer 
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speed  is  the  sum  of  its  eigenvalues.   It  can  be  easily  shown 
that  this  sum  is  numerically  equal  to  the  trace  (the  sum  of 
the  diagonal  elements)  of  the  F  matrix. 

As  our  constraint  on  observer  speed  we  therefore  set  a 
limit,  T,  to  the  most  negative  value  that  we  will  allow  the 
trace  of  the  matrix  £  to  assume.   To  observe  the  effect  of 
T  on  the  resulting  design  we  will  allow  T  to  take  on  the 
values  -9,  -12,  and  -15.   We  see  that  for  our  initial  esti- 
mate for  G  that  the  trace  of  F  is  -9  and  that  this  G  will 
therefore  satisfy  any  of  the  values  of  T  previously  stated. 

The  design  algorithm  of  Section  IID  was  used  to  compute 

* 

the  optimal  observer  gain  matrix  G  ,  the  corresponding  value 

of  the  performance  measure  J(G  ),  the  observer  eigenvalues 
X-.    and  X„,  and  the  worst  case  initial  condition  error  vector 
for  various  choices  of  the  design  specifications  Q^   and  r  in 
the  performance  measure  (20)  and  for  various  choices  of  the 
constraint  T.   The  plant  and  observer  were  then  simulated 
on  the  IBM  360  digital  computer  for  each  design  specifica- 
tion  using  the  corresponding  G   and  the  worst  case  initial 
condition  error. 

Table  I  shows  the  effect  of  varying  the  design  parameter 

*     * 

r  on  G  ,  J (G  ),  and  the  observer  eigenvalues  A^  and  X-.   For 

these  cases  the  matrix  ^   is  fixed  at 


a  = 


1   0 
0   1 


and  T  is  fixed  at  -12.   Note  that  in  each  case,  the  value 
for  g-  is  zero.   Since  the  trace  F_   must  be  greater  than  or 
equal  to  T  we  see  that 
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-  12  -  g^  >  -  12  ; 

is  satisfied  only  for  g^  less  than  or  equal  to  zero.  The 
fact  that  g,  is  always  zero  indicates  that  the  constraint 
is  alvrays  active. 

The  effect  of  varying  the  parameter  r  on  the  observer 
error  time  response  is  shown  in  Fig.  3.   Note  that  as  the 
value  of  r  is  increased,  the  estimate  error  for  both  states 
reaches  zero  sooner.   For  the  estimate  error  for  the  state 
x,(t),  however,  this  is  at  the  expense  of  larger  error 
values  in  the  early  part  of  the  operating  interval. 

In  Table  II  the  effect  of  the  design  parameter  q..^  on 

*     * 

G  ,  J (G  ),  X^  and  A   is  shown.   Here,  the  value  of  r  is 

fixed  at  zero  and  the  value  of  T  is  fixed  at  -12.   The  ma- 
trix 1^  has  the  form 


a  = 


^11 


0" 


Again  we  note  that  the  constraint  T  is  active. 

The  observer  error  time  response  is  shown  in  Fig.  4. 
The  effect  of  increasing  q...  on  the  plant  states  is  similar 
to  that  of  increasing  r,  although  somewhat  less  dramatic. 

In  Table  III  the  effect  of  the  design  parameter  q^^  on 
G  ,  J  (G  ) ,  X^,   and  A_  is  shown.   Again,  the  value  of  r  is 
fixed  at  zero  and  the  value  of  T  is  fixed  at  -12.   The  ma- 
trix 2.  has  the  form 


1 

0 


0 


'22 
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«1 

JCG*) 


0 
0.0 

-35.0 

1.04 
-10.9 
-  1.1 


0.0 

19.4 

0.40 

8.9 

3.1 


0.0 

27.5 

0.11 
-6.0-J6.2 
-6.0+J6.2 


Table  I.   The  effect  of  r  on  G  ,  JCG  ),  X   and  X2. 

'1   0' 


a  = 


0  1 


,  T  =  -12 


^11 


Si 

* 

22 

* 

J(G  ) 


0.0 

35.0 

1.04 
10.9 
1.1 


10 


0.0 

30.0 

1.50 
10.4 
1.6 


0.0 

-8.9 

3.  58 
■6.0-J1.4 
•6,0  +  jl.4 


Table  II.   The  effect  of  q.^  on  G  ,  J (G  ),  X^    and  X^ 

r  =  0,  T  =  -12 
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.  KEY: 


T 

1.5       2.0 
TIME(sec) 

(a)   Estimate  Error  for  Plant  State  x-,(t) 


0.0    - 


1.0    - 


CM 

<<dI 


2.0    _ 


-3.0    _ 


-4.0    _ 


TIMECsec) 
1.5  2.0 


(b)   Estimate  Error  for  Plant  State  x,(t) 


Figure  3.   Observer  Estimate  Error  Time  Response  as  a  Function 
of  Time  Multiplier  Exponent  r. 
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1  qii=l 

2  qii=2 

3  q-^;^=10 


f   .      ^   .  TIME (sec) 

UJ   Estimate  Error  for  Plant  State  x^(.t) 


1.5       2.0 

( 
2' 


TIME(sec) 


(b)   Estimate  Error  for  Plant  State  x^(t) 


Figure  4.   Observer  Estimate  Error  Time  Response  as  a  Function 
of  q^^. 
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Once  again  we  see  that  the  constraint  T  is  active. 

The  observer  error  time  response  is  shown  in  Fig.  5. 
The  effect  of  <\yy    on  the  response  is  quite  different  from 
the  effects  of  r  and  q^,  as  we  might  expect.   By  placing 
greater  emphasis  on  the  cost  associated  with  states  which 
represent  higher  derivatives  we  do  not  allow  them  to  devi- 
ate as  far  from  zero  and  thus  tend  to  slow  the  convergence 
process  down.   Hence,  as  q-^  increases,  the  observer  error 
states  both  reach  zero  at  later  values  of  time. 

It  should  be  pointed  out  that  all  observer  configura- 
tions were  simulated  for  a  period  of  two  seconds  for  the 
sake  of  uniformity  in  the  graphical  presentations.   If  the 
response  of  Fig.  5  were  continued  for  a  greater  length  of 
time,  the  curves  would  all  eventually  approach  zero. 

In  Table  IV  the  effect  of  the  design  constraint  T  on 

*     * 

G  ,  J(G  ),  \^    and  A_  is  shown.   Again,  r  is  fixed  at  zero 

and  Q^  is  fixed  at 

"l  o' 


^  = 


0  1 


By  noting  the  values  for  g^  we  see  that  the  constraint  T  is 
always  active. 

The  observer  error  time  response  is  shown  in  Fig.  6. 
We  see  that  the  effect  of  making  T  more  negative  is  to  im- 
prove the  shape  of  the  response  at  the  earlier  values  of 
time.   Note  that  in  the  case  where  T  has  the  value  -9  that 
the  observer  error  in  state  x^(t)  actually  gets  worse  for 
early  values  of  time. 
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^2  2 

1 

2 

10 

* 
Si 

0.0 

0.0 

0.0 

* 
82 

-35.0 

-38.5 

-43.2 

J(G*) 

1.04 

1.46 

3.21 

^1 

-10.9 

-11.2 

-11.7 

^2 

-    1.1 

-    0.8 

-    0.3 

Table  III.   The  effect  of  o^^^   on  G  ,  J(G  ),  X^  and  ^2 
r  =  0,  T  =  -12 


Si 

* 

S2 

J(G*) 


-9 


-3.0 

27.1 

8.30 
•4.5-J4.2 
■4.5  +  J4/2 


12 


0.0 

-35.0 

1.04 
-10.9 
-  1.1 


15 


3.0 

-44.9 

0.16 
-11.8 
-  3.2 


Table    IV.      The    effect    of  T   on   G    ,    J (G   ),    X      and   X 
fl    0" 


0    1 


,    r   =    0 
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<a)| 


1.5        2.0 

TIMECsec) 

-(a)   Estimate  Error  for  Plant  State  x„ (t) 


TIME(sec) 
1.5       2.0 

L 


(b)   Estimate  Error  for  Plant  State  x_(t) 


Figure  5.   Observer  Estimate  Error  Time  Response  as  a  Function 
°f  ^22- 
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4J 


KEY 

1  T=-9 

2  T=-12 

3  T=-15 


1.5         2.0 
TIME (sec) 


(a)   Estimate  Error  for  State  X2(t) 


0.0 


<  0)1 


-2.0  - 


-4.0  _ 


TIME  (sec) 
1.5       2.0 


(b)   Estimate  Error  for  State  y^T^it) 


Figure  6.   Observer  Estimate  Error  Time  Response  as  a  Function 
of  Constraint  T. 
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IV.   OBSERVER/CONTROLLER  DESIGN 

A.  TRADITIONAL  DESIGN  TECHNIQUE 

The  traditional  method  for  the  design  of  feedback  con- 
trol systems  utilizing  observer  estimates  of  unmeasured 
states  is  to  specify  desired  plant  performance  in  the  manner 
described  in  Section  IIA  as  if  all  states  are  available 
through  direct  measurements.   A  feedback  network  or  controller 
is  then  constructed  which  causes  the  closed-loop  system  to 
satisfy  these  performance  specifications.   The  controller 
has  available  as  inputs  either  the  observer  estimate  of  the 
entire  state  vector  or  the  measured  plant  states  and  the  ob- 
server estimate  of  the  remaining  states.   The  observer  is 
designed  to  be  stable  and  to  have  a  characteristic  equation 
with  eigenvalues  which  satisfy  the  somewhat  nebulous  require- 
ment that  the  observer  by  "slightly  faster  than  the  plant" 
[2]. 

In  this  section,  an  alternative  method  is  described  in 
which  the  observer  and  controller  are  designed  simultaneously 
using  a  performance  measure  similar  in  form  to  (2)  in  Sec- 
tion IIB.   This  results  in  a  configuration  in  which  both  the 
controller  design  and  the  observer  design  are  optimal  with 
respect  to  the  specified  performance  measure. 

B.  FULL-ORDER  OBSERVHR/CONTROLLER  DESIGN 

The  full-order  observer  for  the  linear,  time  -  invariant 
system 
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x(t)  =  Ax(t)  +  Bu(t) 

(21) 
)r(t)  =  Cx(t) 

is,  as  shovm  in  Appendix  D, 

zCt)  =  (A-GC)z_(t)  +  G};^Ct)  +  Bu(t)  .     (22) 
We  define  the  control  law  as 

u(t)  =  Kx(t) 

where  K  is  a  constant  Ixn   matrix. 

It  can  be  shown  [2]  that  the  stability  of  a  closed-loop 
system  using  an  observer  output  as  a  substitute  for  the 
plant  states  in  the  control  law  is  unaffected  by  a  stable 
observer.   The  control  can  therefore  be  approximated  by  the 
control  estimate,  u(t) ,  where 

u(t)  =  Kz(t)  .  (23) 

The  structure  for  this  full-order  observer/controller  is 
shown  in  Fig.  7. 

Defining  the  error  state,  £(t)  ,  as 

e(t)  ^  z(t)  -  x(t)  (24) 

and  differentiating  with  respect  to  time  yields 

e(t)  =  z^(t)  -  x(t)  . 

Substituting  (21)  and  (22)  for  x(t)  and  z_(t)  ,  respectively, 
and  solving  for  e(t)  yields 

e(t)  =  (A-GC)e(t).  (25) 

Substituting  (23)  into  (21)  yields 

x(t)  =  Ax(t)  +  BK£(t)  .  (26) 
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Figure  7.   Full-Order  Observer/Controller, 
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Solving  (24)  for  z_Ct)  and  substituting  into  (26)  yields 


x(t)  =  (A+BK)x(t)  +  BKe(t) 


Defining  the  augmented  state  vector,  C(t),  as 


(27) 


?(t)  ^ 


'x(t) 
e(t) 


and  using  (25)  and  (27),  it  is  easily  shown  that 


i(t)  = 


A+BK     BK 
0      A-GC 


C(t) 


(28) 


The  performance  measure  can  be  defined  as 


J  =  r  t^[x'^(t)qx(t)  +  e'^(t)Re(t)]dt. 


(29) 


This  performance  measure  is  chosen  because  it  accommodates 
the  specification  of  convergence  properties  of  both  the 
closed-loop  system  states  and  the  observer  states  and  im- 
plicitly limits  the  expenditure  of  control  effort.   Equa- 
tion (29)  can  be  written  in  terms  of  ?(t)  as 


00 


(t) 


9.9. 

0  R 


C(t)dt. 


(30) 


Equations  (28)  and  (30)  are  in  the  same  form  as  (1)  and 
(2),  respectively.   It  has  been  shown  [2]  that  the  matrix 
4^,  defined  as 

A+BK    BK 


^  ^ 


0 


A-GC 
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has  as  its  eigenvalues,  the  eigenvalues  of  A+BjC  and  the 
eigenvalues  of  A-G£.   The  matrix  H^  is  therefore  a  stability 
matrix  if  both  A+BK  and  A-G£  are  stability  matrices.   It 
can  be  shown  that  the  matrix 

0   R 

is  positive  definite  if  both  2.  a^^d  R  are  positive  definite 
[7] .   Therefore  (30)  can  be  rewritten  as 


J  =  (-l)^''^r:i^(0)V^^^iCO) 


C31) 


where,  by  using  (6) 


[A+BK   BK  ■ 
0     A-GC 


V    +  V 
-s+1   -s+1 


A+BK    BK  ■ 
0     A-GC 


=   Is        s  =  l,2,.  .  .  ,r 


A+BK    BK 
0     A-GC 


-iT 


^1 


^^1 


A+BK    BK 
0     A-GC 


0   R 


Furthermore,  the  algorithm  described  in  Section  IID  can 
be  applied  to  find  K   and  G  ,  the  feedback  and  observer 
gain  matrices  which  are  optimal  with  respect  to  the  per- 
formance measure  (29) . 

A  slight  modification  to  the  design  algorithm  causes  a 
significant  reduction  in  the  number  of  computations  re- 
quired on  each  iteration.   To  see  why,  we  note  that  if  the 
observer  estimate  initial  condition  vector  is  assumed  to  be 
the  null  vector,  that  is 
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z(0)  =  0 


then,  from  (23),  eCO)  becomes 


eCO)  =  -xCO) 


and  the  augmented  state  vector  C(0)  becomes 


i(0)  = 


x(0) 
-x(0) 


(32) 


Partitioning  the  symmetric  matrix  V  _^-.  as 


-r+1 


^11  -Lii 


T 
V    V 
-12   -2  2 


where  each  element  is  an  nxn  matrix,  then  substituting  (32) 

and  (33)  into  (31)  yields 

(34) 

J  =  (-l)'"''^r:J(0)[V^^-V^2-Xi2^^V22]x(0). 

We  now  define  the  nxn  matrix  V.^  as 

—r+1 


V;,i  =  [Vii-V,2-V,2^-V22] 


(35) 


Therefore 


J  =  (-1)^^^:  J(0)v;^^x(0) 


(36) 


The  matrix  V'  .,  has  n  eigenvalues,  while  the  matrix  V  ^, 
—r+1         ^        '  —r+1 

has  2n  eigenvalues.   We  see,  therefore,  that  the  effort 
required  in  finding  J,  can  be  significantly  reduced  by  using 
(35)  and  (36)  in  step  three  of  the  design  procedure. 
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C.   REDUCED-ORDER  OBSERVER/CONTROLLER  DESIGN 
For  the  linear,  time -invariant  system 

x(t)  =  Ax(t)  +  BuCt) 
xCt)  =  Cx(t), 

we  define  the  control  law  _,  _ 

u(t)  ^  KxCt) 


(37) 


(38) 


where  K  is  a  constant  Zxn   matrix.   We  define  the  vector 


x(t)  ^ 


rcty 

w(t) 


(39) 


where  w(t)  is  an  (n-m)xl  vector.   We  let 

w(t)  =  Dx(t) 

where  B_   is  an  (n-m)xn  matrix,  selected  so  that 

_  _  1 
C 

D 

exists.   By  proper  selection  of  the  state  variables,  C^  will 
have  the  form 

C  =  [  iiO  ]  . 

A  notationally  convenient  choice  for  D  is 

D  =  [  0  :;^  ]  . 

For  this  choice  of  C  and  D,  w(t) ,  the  estimate  for  w(t) , 
can  be  obtained  by  constructing  a  reduced-order  observer, 
as  developed  in  Appendix  D,  whose  state  equations  are  given 
by 
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zCt)  =  (A22-GA^2^^'^^^  ""  <^A2l"^ll'''^^22"^12^^^^^^^ 


+  (B2-GB^)u(t) 


(40) 


where 


w(t)  =  zCt)  +  G^Ct). 


(41) 


The  matrix  K  and  the  vector  x(t)  can  be  partitioned  as 


K  =  [K^:K°]  ,   x(t)  = 


where,  for  the  stated  choice  for  C  and  D, 


x'^Ct) 
i°Ct) 


)r(t)  =  x^'Ct),  w(t)  =  x°Ct),  wCt)  =  £°(t). 

The  control  law,  given  by  (38),  can  therefore  be  approximated 
by 


u(t)  =  K 


x"^(t) 
i°(t) 


or 


G(t)  =   [k'":K°] 


iCt) 
z(t)  +  Gx(t) 


or 


u(t)  =  K\(t)  +  K°z  (t)  +  K°Gx(t)  . 


(42) 


The  structure  for  this  reduced-order  observer/controller  is 
shown  in  Fig.  8. 

We  define  the  observation  error,  e_(t)  ,  as 

e(t)  ^  w(t)  -  w(t).  (43) 

Substituting  (41)  for  w(t)  and  (39)  for  w(t)  yields 


e_(t)  =  _z(t)  +  Gjr^Ct)  -  Dx(t). 


(44) 


40 


xCt) 


D=*   J 


7* 


B 


xCt) 


zct) 


p 


^2r^ll"fA22-GA^2^G 


^ 


^=- 


/ 


=> 


B^-GB^ 


z(tl 


+ 


u(t) 


N  / 


zCt) 


^22-^12 


<± 


K° 


<^ 


^m 


C^ 


"^^ 


Figure  8.   Reduced-Order  Obscrver/Controll 


er, 
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Substituting  (37)  for  };;(t)  and  rearranging  yields 

e(t)  =  z(t)  +  (GC-D)x(t)  . 
Differentiating  with  respect  to  time  yields 

e(t)  =  iCt)  +  (GC-D)x(t)  . 

Substituting  (37)  and  (40)  for  x(t)  and  z_(t)  ,  respectively, 
and  solving  for  £(t)  yields,  as  shown  in  Section  IIIB, 

|_(t)  =  (A22-GA^2^^'^^^' 
Solving  (44)  for  z^(t)  +  G^(t)  and  substituting  into  (42) 
yields 

u(t)  =  K"^x(t)  +  K°(e(t)  +  Dx(t)). 

We  note,  however,  that 

Dx(t)  =  x°(t) 

and 


m, 


^(t)  =  x"'(t) 


Therefore , 


u(t)  =  Kx(t)  +  K"e(t) 


(46) 


Substituting  this  estimate  for  u(t)  into  (37)  yields 


x(t)  =  Ax(t)  +  B(Kx(t)  +  K''e(t) 


which,  after  rearranging,  becomes 


x(t)  =  (A+BK)x(t)  +  BK"e(t) ) 


(47) 


We  define  the  augmented  state  vector,  C(t),  as 


i(t) 


x(t) 
e(tj 
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Using  this  definition  and  the  expressions  for  x^t)    and  _e(t) 
given  by  (47)  and  (45),  respectively,  yields 

o 


Br 


,  i(t)  = 


A+BK 
0 
The  performance  measure  can  be  defined  as 


-22"— 12 


e(t)  . 


(48) 


(49) 


=  J    t^[x'^(t)2x(t)  +  J(t)Re(t)]dt 


which  can  be  rewritten  in  terms  of  the  augmented  state  vec- 
tor, K_(t)  ,  as 


=  /  tV(t) 


0   R 


^(t)dt.       (50) 


Equations  (48)  and  (50)  are  in  the  same  form  as  (1)  and 
(2),  respectively.   Furthermore,  the  matrix 

A+BK    BK° 

i    ^22-^12 

is  a  stability  matrix  if  both  A+BK  and  A„^-GA^~  are  stability 
matrices.   Also,  the  matrix 


0   R 


is  positive  definite  if  both  Q^   and  R  are  positive  definite. 
The  performance  measure  (50)  can  therefore  be  rewritten 
as 


J  =  (-l)''^^r:i'^"(0)V^^jC(0) 


(51) 


4  3 


where,  by  using  (6), 


A+BK   BK' 


Li  hl-^llA 


^s.l 


^s.l 


A+BK    BK 
0 


o 


-22"— 12. 


V 
— s 


s  =  1, 2  , . . .  ,r 


A+BK    Br 


L-0   hl-^ll. 


^1 


^^1 


A+BK 


0 


BK' 


-2  2'— 12. 


0  R 


Furthermore,  the  design  procedure  of  Section  IID  can  be 

*      * 

applied  to  determine  K   and  G  ,  the  feedback  and  observer 

gain  matrices  which  are  optimal  with  respect  to  the  speci- 
fied performance  measure  (49) , 

Once  again,  it  is  possible  to  reduce  the  number  of 
computations  required  for  each  iteration  of  the  design 
algorithm.   In  this  case,  we  set  the  observer  estimate, 
w(0)  =  £,  which  yields  e(0)  =  -Dx(0).   Therefore  _§^(0)  be- 
comes 

^  x(Or 
■DxCO) 


ao)  = 


Substituting  this  into  (51)  yields 


(52) 


(-1)^^^: 


x(0) 
-Dx(0) 


Vl 


x(0)  - 
-Dx(0) 


Partitioning  the  (2n-m) x (2n-m)  symmetric  matrix  V  ^.  as 


-r+l 


^11   ^12 

T 

V     V 
-12   -2  2 
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where  V^ ,  is  an  nxn  matrix,  V^ 2  is  an  nx(n-m)  matrix,  and 

V~2  is  an  Cn-m)x(n-m)  matrix  and  substituting  into  (52) 

yields,  after  expanding 

(53) 

J=  (-l)^''-^r:x'^(0)  [V^-^-d'^V^2'^-^12°^^'^^22-^-'^°^' 

We  define  the  nxn  matrix  V'  ^  as 

— r+1 


V;,i  ^  ilil-^^U      -   ^12°  '   ^''^22^ 


(54) 


Therefore,  (53)  can  be  rewritten  as 
J  =  (-l)'"''^r:x^(0)V;^^x(0) 


(55) 


The  matrix  V'  ,  has  n  eigenvalues,  while  the  matrix  V  ,, 
—r+1         ^        '  —r+1 

has  2n-m  eigenvalues.   We  see,  therefore  that  the  effort 
required  in  determining  J.  in  step  three  of  the  design  pro- 
cedure can  be  reduced  by  using  (54)  and  (55). 

D.   EXAMPLE  OBSERVER/CONTROLLER  DESIGN  PROBLEM 

To  illustrate  the  application  of  the  design  algorithm 
to  observer/controller  design  we  consider  the  plant  whose 
transfer  function  is  given  by 

Y(S) 


U(S)    S2  +  15S  +  50  • 

The  state  and  output  equations  for  this  plant  can  be  written 
as 


x(t)  = 


0    1 

■50  -15 


x(t) 


u(t) 


i:(t)  =  [10  ]x(t) 


Here , 
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Then 


A  = 


0 

1 

'    B    = 

0 

50 

-15 

1 

,  C  =  [  1  0  ] 


A+BK  = 


0 

1  ' 

+ 

0 

50 

-15 

_1_ 

[k,  k^] 


|kj^-50   k2-15 


Also 


where 


F   A22  -  GA^2 


A22  =  [-15],   A 


12 


[1]  ,  G  =  [g] 


Therefore 


F  =  -[15+g]. 


We  must  select  an  initial  estimate  for  both  K  and  £ 
which  causes  the  matrices  A+M  and  IF  to  be  stability  matrices 
As  in  the  case  of  the  observer  alone  we  must  constrain  sys- 
tem speed.   In  this  case,  however,  we  sum  the  eigenvalues 
of  the  closed  loop  system.   It  is  easily  shown  that  this 
sum  is  given  by  the  trace  of  the  matrix 

,0 


A+BK 


0 


BK 


^22'— 12 


which  is  given  by 


k2  -  30  -  g 
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We  arbitrarily  select  the  trace  constraint  T  to  have  a  value 
of  -25.   Therefore 

-25  <  k2  -  30  -  g. 

We  see  that  the  open-loop  system  is  stable.   We  therefore 
choose  an  initial  estimate  for  K  of  [0  0] .   We  must  now 
choose  a  value  for  g  which  satisfies  both 

-25  <  -  30  -  g 
and 

-(15+g)  <  0  . 

Or,  combining  the  two  inequalities, 

-15  <  g  <  -  5. 

We  choose  for  our  estimate  of  G,  therefore, 

G  -  [-10]. 

The  design  algorithm  of  Section  IID  was  used  to  compute 
the  optimal  controller  feedback  gains  K   and  observer  gains 
G  ,  the  corresponding  value  for  the  performance  measure 
J  (K  ,G^  ),  the  plant  and  observer  eigenvalues  X^,  X„  and  X_, 
and  the  worst  case  plant  and  observer  estimate  error  ini- 
tial conditions  for  various  choices  of  ^  and  r  in  the  per- 
formance measure  (49).   The  constraint  T  was  fixed  at  -25 
and  the  design  parameter  R  was  fixed  at 

R  =  [1] 

in  all  cases.   The  plant  and  observer  were  simulated  on  the 
IBM  360  digital  computer  for  each  performance  measure  spe- 
cificd  using  K^  ,  G    ,    and  the  corresponding  plant  and  esti- 
mate worst  case  initial  conditions. 
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Table  V  shows  the  effect  of  varying  the  design  parameter 

*    *     *   * 

r  on  K  ,  G^  ,  J(K  ,G  )  ,  A^  ,  X-,  and  X_.   In  this  case  the  Q^ 

matrix  is  also  fixed  at 


9.  = 


1      0 
0   1 


Figure  9  shows  the  effect  of  varying  the  design  parameter  r 
on  the  closed-loop  system  response  to  the  worst  case  initial 
conditions.   We  see  that  increasing  r  improves  the  system 
settling  time  in  both  states  at  the  expense  of  larger  values 
of  x-(t)  in  the  early  part  of  the  operating  interval. 

Table  VI  shows  the  effect  of  varying  the  design  parameter 
q..  on  K  ,  £  ,  J (K  ,G  ) ,  X- ,  A„,  and  X_.   In  this  case  r  is 
fixed  at  zero  and  Q^  has  the  form 


S  = 


■qil   0 
0    1 


Figure  10  shows  the  effect  of  varying  q...  on  the  closed-loop 
system  response  to  the  worst  case  initial  conditions.   We 
see  that  as  the  value  of  q.,^  is  increased,  the  settling  time 
improves  for  both  states  at  the  expense  of  larger  values  of 
x_(t)  in  the  early  part  of  the  operating  interval. 

Table  VII  shows  the  effect  of  varying  the  design  parame- 
ter  q^^  on  K  ,  G  ,  J (K  ,G  ),  \^,    X^,    and  A^.   In  this  case 
r  is  fixed  at  zero  and  ^   has  the  form 

1    0 

0   ^l22 
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k 
k 
g 


J(K    ,G    ) 


0 

25.7 

-  9.4 

-14.4 
1.02 

-  0.6 
-23.3 

-  1.0 


17,1 
9.0 

14.0 
0.34 
1.0 

20.8 
3.2 


-217.4 

-  8.0 

-13.0 
0.06 

-  2.0 
11.5-J11.6 
11.5+J11.6 


Table   V.      The    effect   of   r   on   K* ,    G * ,    J(K*,G*),    X^,    X^,    and   X^. 
-^    °       ,    T    =    -25,    R   =    [1]  . 


a  = 


0    1 


^11 

1 

2 

10 

* 

^1 

25.7 

-15.5 

-25.9 

* 

-    9.4 

-    9.5 

-    9.0 

is 

g 

-14.4 

-14.5 

-14.0 

J(K    ,G    ) 

1.02 

1.46 

3,37 

h 

-    0.6 

-    0.5 

-    1.0 

h 

-23.3 

-23.0 

-20,3 

h 

-    1.0 

-    1.5 

-    3.7 

Table  VI.   The  effect  of  q^^  on  K  ,  g",  J (K  ,G  ),  X^,    X^, 


and  X 


0,  ^T 


25,  R  =  [\] 
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KEY 

1  r=0 

2  r=l 

3  r=2 


TIME(sec) 
(a)   Closed-Loop  System  State  x^(t) 


TIME(sec) 
2.0 


(b)   Closed-Loop  System  State  x^ (t) 


Figure  9.   Closed -Loop  System  Time  Response  as  a  Function  of 
Time  Multiplier  1-xponcnt  r. 
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KEY 

i.qn=i 

2  q^^=2 

3  qi;L"^° 


X 


1.5        2.0 
TIME (sec) 
(a)   Closed-Loop  System  State  x^  (t) . 


0.0 


1.0- 


X 


2.0- 


0.5 


1.0 


TIME(sec) 
1.5         2.0 


(b)   Closed-Loop  System  State  x^(t) 


Figure  10. 


Closed-Loop  System  Time  Response  as  a  Function 
of  q^^. 
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''22 

1 

2 

10 

< 

25.7 

32.7 

42.2 

* 

-    9.4 

-    9.4 

-    9.6 

* 

s 

-14.4 

-14.4 

-14.6 

*      * 

J(K    ,G    ) 

1.02 

1.44 

3.18 

h 

-    0.6 

-    0.6 

-    0.4 

h 

-23.3 

-23.7 

-24.3 

h 

-    1.0 

-    0.7 

-    0.3 

*       * 


Table   VII.       The    effect    of   q^^    °^   ^    '£   '^C^    >1  )>    -^i*    ^7 


and   X 


r=0, 


T=-25,    R=[lT 
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Figure  11  shows  the  effect  of  varying  q-o  on  the  closed-loop 
system  response  to  the  worst  case  initial  conditions.   We 
see  that  as  q^^  is  increased,  the  system  becomes  more  slug- 
gish, requiring  longer  settling  times,  as  we  would  expect. 


53 


KEY 


1.0 


.0.5  _ 


r- 

X 


0.0 


TIME  (sec) 
(a)   Closed-Loop  System  State  x^ (t) 


0.0 


-0.5  _ 


TIMECsec) 


Closed-Loop  System  State  Xp(t) 


Figure  11.   Closed-Loop  System  Time  Response  as  a  Function 


of  q 


22 
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V.   CONCLUSIONS 

A  new  technique  for  observer  and  observer/controller 
design  which  requires  the  minimization  o£  an  integral  of 
time-multiplied-square -error  (ITSE)  performance  measure  has 
been  developed.   An  algorithm  which  is  well  suited  to  digi- 
tal computation  and  which  is  especially  useful  for  higher- 
order  systems  has  been  presented.   An  example  of  a  computer 
program  which  implements  the  algorithm  for  the  case  of  an 
observer/controller  is  shown  in  Appendix  E.   The  computa- 
tions are  performed  off-line  with  a  typical  computation  time 
of  thirty  seconds. 

It  has  been  demonstrated  through  the  use  of  examples 
that  modification  of  the  time  response  of  an  observer  or 
a  closed-loop  system  utilizing  an  observer/controller  to 
initial  condition  error  can  be  achieved  by  a  suitable  change 
in  the  appropriate  design  parameters  of  the  ITSE  performance 
measure.   Furthermore,  it  has  been  demonstrated  that  the 
qualitative  change  in  the  system  time  response  reflects 
directly  the  qualitative  change  in  the  performance  measure 
design  parameter. 

It  is  recognized  tliat  for  second-order  and,  perhaps, 
third-order  s ingle  -  input ,  single-output  systems  the  rela- 
tionships between  traditional  design  specifications  and 
system  transfer  functions  are  simple  and  relatively  easy  to 
use.   In  these  cases,  the  design  procedure  suggested  above 
is  generally  not  efficient.   However,  when  dealing  with 

5  5 


higher-order,  multiple-input,  multiple-output  systems,  the  ad- 
vantages o£  the  new  technique  become  more  significant.   For 
these  systems,  the  design  specification  and  system  transfer 
function  relationships  are  highly  complex  and  the  traditional 
methods  are  therefore  extremely  difficult  to  use.   However, 
the  new  technique  offers  a  relatively  straightforward  ap- 
proach to  design  vvhich  does  not  become  more  complex  as  the 
system  order,  the  number  of  inputs,  or  the  number  of  outputs 
increases. 
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APPENDIX  A 

00 

EVALUATION  OF  THE  INFINITE  TIME  INTEGRAL  /  t^£(t)gs (t) dt 

o 

For  the  system  of  linear,  first-order  differential 
equations 

s(t)  =  HsCt)  (lA) 

if  there  exists  a  scalar  function 

V^(s(t))  =  /(t)V^s(t)  C2A) 

subject  to  the  conditions 

1)  V^(0)  =  0 

2)  V^(sCt))  >  0,  for  all  s_(t)  ^  0_ 

3)  V^(s_(t))-  00,  as  I  |^(t)|  I  ->  oo 

4)  V^(s(t))  <  0,  for  all  s^(t)  7^  0 

5)  V^(s(t))  =  0,  £(t)  -  0 

then  the  system  (lA)  is  asymptotically  stable  [4] . 

Since  V,(s_(t))  is  a  scalar,  it  is  equal  to  its  transpose 
Hence,  from  (2A) ,  we  see  that  the  matrix  V^  is  symmetric. 
We  now  consider 

V^(s_(t))  =  i'^Ct)V^sCt)  +  s'^(t)V^i(t). 
Substituting  (lA)  for  s_(t)  yields 

V^(s(t))  =  s'^(t)H^V^s(t)  +  s'^(t)VjHs(t) 

which,  after  rearranging  terms  becomes 
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V^Cs(t))  =  s'^(t)[H'^V^  +  V^H]s.(t). 
We  now  define  the  matrix  ^   as 

and  note  that  since  V,  is  symmetric,  then  ^  must  also  be 
symmetric.   Therefore 

V^CsCt))  =  -Jct)^(t).  (3A) 

At  this  point  we  note  that  equation  (2A)  insures  that 
conditions  one  and  three  are  satisfied.   Condition  two  is 
satisfied  if  V^  is  positive  definite  and  conditions  four 
and  five  are  satisfied  if  ^  is  positive  definite. 

Integrating  the  left  hand  side  of  (3A)  over  the  finite 
time  interval,  T,  yields 

T  T 

rv^(s(t))dt  =  r  dv^cs(t)) 

o  o 

T 

Jv^(sCt))dt  =  V^(s(T))  -  V^(s(0)). 
o 

If  (lA)  is  asymptotically  stable,  then 
lim  V-^(£(T))  =  0. 


Therefore , 


00 

j'v^(s_(t))dt    =    -V-^(s(0)) 


Or,  from  (3A) , 


j   s'^(t)gs_(t)dt  =  s^(0)V^s_(0) 


(4A) 


where  V^  is  the  solution  to  the  algebraic  matrix  equation 


H  li  +  V^H  =  -^ 


(5A) 


We  now  consider 


^[ts(t)V^s(t)]  =  s'^(t)V^s(t)  +  t[s'^(t)V^s,(t)  +  J(t)V^sCt)] 


V^(sCt))  +  tV^(s(t)) 


Hence 


CO  00  00 

/  dT^^^^'^^^^l^'^^^^  ^  /  V-LClCt))dt  +  j    tV^(s(t))dt 
o  o  o 

It  can  be  shown  that  if  s_(t)  is  the  state  vector  of  an 

asymptotically  stable,  free  system,  then  for  a  finite  n. 


d[t"s_'^(t)V^s_(t)]    =    t^s'^(t)V^£(t) 


=    0 


Therefore 


0    = 


J      V^(£(t))dt    +    J      tV^(^(t))dt 
o  o 


or 


CO  00 

J    tv^(srt))dt  =  -J    Vj(^(t))dt 
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or,  substituting  (3A)  for  Vj^(s_Ct)), 


oo  oo 

J    ts'^(t)qs_Ct)dt  =  -  J   V^(s(t))dt. 


(6A) 
o 


Now 


suppose  there  exists  a  scalar  V^(s_(t))  such  that 

.  V2(s(t))  ^  s^(t)V2s(t)  C7A) 

and 

V2Cs(t))  =  V^(s(t)).  C8A) 


Then 


or 


V2Cs(t))  =  i'^(t)V2s(t)  +  s'^(t)V2i(t) 


V2Cs(t)3  =  s'^(t)[H'^V2  +  V2H]s(t) 


However,  from  (8A)  and  (2A)  we  have 


V2(s(t))  =  s'^(t)V^s(t) 


Therefore,  we  see  that 


h'^V2  ^   V2H  =  V^ 


where,  since  V..  is  symmetric,  then  V^  is  symmetric.   From 
(8A)  we  have 


J  V^Cs.(t))dt  =  -  j  V2(l(t))dt 
o  o 


00 

V^(s(t))dt  =  -  j    dV2(s^(t)) 
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-J  V^Cs(t))dt  =  V^CsCO)) 
o 

00 

-  r  V^Cs(t))dt  =  J(0)V2S(0) 


Substituting  (6A)  yields 


CO 

-f  ts^^Ct)^(t)dt  =  s'^(0)V2s(0) 


(9A) 


wh 


ere  V^  is  the  solution  to  the  algebraic  matrix  equations 


h'^V^  -   V^H   =  V^ 


h'^v^  +  v^H  =  -  g. 


(lOA) 


The  argument  can  be  continued  for  higher  powers  o£  t. 
Thus  for  a  stable  system  (lA)  and  a  symmetric,  positive 
definite  ^, 

00 

j    tVCt)gi(t)dt  =  (-l)^^lr:s^(0)V^^^s(0)   (llA) 
o 

where  V  ^ ,  is  the  solution  to  the  simultaneous  algebraic 
— r+1  ^ 

matrix  equations 

h'^V  ^i+V^,H  =  V        s  =  l,2,...,r 
—  — s+1    — s+1—   — s  '  '    ' 


(12A) 


h'^V^  +  V-^H 


61 


APPENDIX  B 

A  METHOD  FOR  THE  DIRECT  SOLUTION  OF  THE  SYMMETRIC  MATRIX  V   , 

— r+1 

A  procedure  for  generating  V   ..  from  the  system  of 
simultaneous  linear  algebraic  matrix  equations 


H  V  ^1  +  V  ^,H  =  V 
s  +  1    — s  +  1—   —J 


J-}^}«  •  •  y    L 


H'V^  +  V^H    =  -  g 

is  described  in  Reference  3.   The  procedure,  which  is  es- 
sentially repeated  here  with  minor  amplification  and  nota- 

tional  changes  involves  the  direct  solution  for  V  ,  -,  and 
'^  —r+1 

bypasses  the  need  for  solving  the  (r+l) (n+l)n/2  simultaneous 
equations. 

We  first  form  the  matrix  S_  as  follows: 

1)  Form  the  symmetric  array  of  integers 

12  4  7 

2  3  5  8 

4  5  6  9 

7  8  9  10 

•n(n+l)/2 

where  n  is  the  dimension  of  H,  Q,  and  V   ,. 

— '  -^ '     —r+1 

2)  Form  the  n  intermediate  square  matrices,  S., 
i=l,2,...,n  whose  elements  arc  equal  to  one  and  whose  order 
is  n(n+l)/2. 


ss 
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3)  Proceed  along  the  first  row  o£  I    and  for  each 
integer,  j,  that  is  missing  from  the  sequence 

{1,2 , . . . ,n(n+l)/2} ,  set  all  entries  in  the  j    row  and  column 
of  S^-j  equal  to  zero. 

4)  Proceed  row  by  row  through  the  H  and  S^  matrices 
simultaneously.   For  the  first  non-zero  entry  in  the  first 
row  of  S^^  ,  enter  the  value  h-^,  for  the  second  non-zero 
entry  in  the  first  row  of  S^.  ,  enter  h..  ^ ,  and  so  forth  until 
all  entries  of  H  have  been  used  to  replace  all  the  non-zero 
entries  of  S^^  . 

5)  Multiply  the  first  column  of  S^.,  by  two. 

6)  Proceed  along  the  second  row  of  I^   as  in  step  three, 
making  the  appropriate  entries  in  the  S^  matrix. 

7)  Repeat  step  four  for  the  £^  matrix. 

8)  Repeat  step  five  for  the  second  column  of  the  S„ 
matrix. 

9)  Continue  the  procedure  until  all  of  the  S*.  matrices 
have  been  generated. 

10)  Form  the  S  matrix  from 


n    T 

s  =  .z,  s.   . 

—    1=1—1 

We  next  form  the  vector  v  .  .,  from 

— r+1 

— r+l    —      -»- 

where  v   -,  and  c^  are  the  leading  diagonal  elements  of  the 
symmetric  v  ^^  and  (^   matrices,  respectively.   That  is. 


6  3 


3^r+l 


'^nl 


V 


12 


V 


22 


V 


13 


V 


23 


33 


V 


nn 


^ 


^11 
^12 
^22 
^13 
^23 
^33 

q 


The  matrix  V  ^i  can  now  be  generated  from  the  vector  v  _,. 
— r+1  ^  — r+1 

in  a  simple  manner. 
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APPENDIX  C 

EVALUATION  OF  THE  QUADRATIC  FORM  s^^  (0)  (- 1)  ^'^■'■r !  V  +;,^s(0) 

Let  X  be  an  n  dimensional,  unit  length  vector.   That  is 

x^x  =  1.  (IC) 

The  n  dimensional  vector  y.  whose  direction  is  the  same  as 
that  o£  X  but  whose  length,  p,  is  arbitrary  is  therefore 
given  by 

I  =  px  •  C2C) 

We  consider  the  quadratic  form 

C  =  x.'^A^  C3C) 

where  A  is  symmetric.   Substituting  (2C)  for  y_   yields 

C  =  p^x'^Ax  .  C4C) 


Let  us  constrain  ^  to  lie  on  or  inside  an  n  dimensional 
hypersphere  of  radius  R.   That  is, 

0  <  p  <  R.  (5C) 

We  assume  we  have  no  further  a   priori    knowledge  of  the 
vector  Y.   ^i""^  therefore  wish  to  consider  its  worst  case  value, 
that  is  the  value  of  y_   which  causes  C  to  take  on  its  maximum 
value  for  a  given  A. 

We  see  from  (4C)  that  C  is  maximized  for  a  given  A  wlien 

T 
p  is  a  maximum  and  when  tlic  scalar  quantity  x  Ax  assumes  its 

maximum  value.   From  (_SC)    wc  see  that  the  maximum  value  for 

p  is  R.   Therefore, 
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C  =   R^   max{x'^Ax}.  (60) 

max  —  —  ^      ^ 

We   define   the   quantity,    q   as 

q   =   niax{x^Ax}  C7C) 

It  can  be  shown  [8]  that  the  matrix  A  in  the  form 

x'^Ax  C8C) 

T  T 

where  x  x  =  1  and  A   =  A,  has  n  real  eigenvalues,  X.  and  n 

corresponding  real  eigenvectors  x-  j  i  -  l,2,...,n.   Further- 
more, the  maximum  value  of  the  form  (8C)  is  the  largest 
eigenvalue  of  A.   That  is 

q  =  maxIx'^Ax}  =  X    [A].  (9C) 

^       —  —     max  "•— ■"  ^      -^ 

Also,  the  vector  which  gives  this  maxim.um  value  is  the  eigen- 
vector corresponding  to  the  largest  eigenvalue. 
Substituting  (9C)  into  (6C)  yields 

C    =  r2  X    [A]  .  (IOC) 

max       max"-— ■•  ^    ^ 

Thus,  C     is  the  value  for  C  in  (3C)  which  represents 
'   max  ^      ^  ^ 

the  worst  case  value  for  ^  for  a  given  A.   We  wish  to  mini- 
mize C  as  a  function  of  A  for  this  worst  case  condition. 
Since  R  is  a  fixed  quantity,  it  plays  no  role  in  the  mini- 
mization of  C  and  we  can  therefore  assume  that  it  is  equal 
to  one.   This  is  equivalent  to  assuming  that  ^  lies  on  the 
unit  hypcrsphere. 

We  are  concerned  with  minimizing  tlie  quantity 


J  =  J(0)f-l)^^^r:V^.^-,sfO)  (IIC) 
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r+1 
where  the  matrix  (-1)    rIV  _^.    is  symmetric  and  positive 

definite  and  where  s_(0)  represents  the  worst  case  initial 

condition  o£  the  state  vector  s_(t)  .   We  can  assume,  without 

loss  o£  generality,  that  s^(0]  lies  on  the  unit  hypersphere. 

From  the  argument  above,  we  see  that  J.,  the  maximum  value 

for  J  for  a  given  parameter  set  {P},  is  given  by 

J,  =  X    [(-Ij^^^rlV  ^-,]  (12C) 

1    max"-  ^  ^  — r+l-"  ^   -* 


and  that  the  worst  case  initial  condition  vector,  s_(0)  ,  is 

r+1 

the  eigenvector  which  corresponds  to  X         [(-1)    rIV  _,-,]. 
&  f  max"^  ^  ^  — r+l-" 
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APPENDIX  D 
OBSERVER  THEORY 

A.   THE  NEED  FOR  OBSERVERS 

Given  a  free  plant  characterized  by  the  state  and  output 
equations 

xCt)  =  Ax(t)  (ID) 

^(t)  =  CxCt)  (2D) 

where  x(t)  is  the  nxl  state  vector,  ^(t)  is  the  mxl  measure- 
ment vector,  A  is  a  known  nxn  matrix,  and  C  is  a  known  mxn 
matrix,  we  wish  to  know  the  value  o£  the  state  xC^) •   If  C 
were  an  nxn  matrix  representing  n  linearly  independent  mea- 
surements, then  C        would  exist  and  x(t)  could  be  found  from 

x(t)  =  C'^i^Ct).  (3D) 

This  is  not  usually  the  case,  however,  so  we  will  assume 
that  there  are  fewer  linearly  independent  measurements  than 
there  are  states.   That  is,  we  will  assume  that  m  is  less 
than  n. 

Since  C^   does  not  exist  for  m  less  than  n,  an  alterna- 
tive scheme  must  be  found  to  determine  x(t) .   One  method  of 
estimating  the  state  vector  is  to  construct  a  model  of  the 
plant  and  operate  it  under  the  same  conditions  as  the  plant. 
Such  a  model  is  called  an  observer.   The  following  sections 
develop  the  state  equations  for  two  different  types  of 
observers . 
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B.   FULL -ORDER  OBSERVERS 

Let  the  observer  state  equations  be  given  by 

_z(t)  =  AzCt).  (4D) 

The  solution  to  (ID)  is  given  by 

x(t)  =  e-^x(O)  (5D) 

and  the  solution  to  (4D)  is  given  by 

^(t)  =  £-^^(0)  .  (6D) 

Therefore,  if  z_CO)  =  xCO)  ,  the  observer  state  will  be  iden- 
tical to  the  plant  state  for  all  t  greater  than  or  equal  to 
zero.  Since,  however,  in  general  z_(0)  does  not  equal  x(0)  , 
the  identical  dynamics  of  (ID)  and  (4D)  would  result  in  the 
persistence  of  the  initial  condition  error  and  the  observer 
state  would  not  provide  an  accurate  estimate  of  the  plant 
state. 

An  alternative  observer  is  one  which  contains  not  only 
the  plant  model  but  also  a  feedback  forcing  function  made 
up  from  an  error  signal.   Since  the  only  accurate  indica- 
tion we  have  of  the  states  is  the  measurement  vector  y;(t)  , 
let  us  assume  that  the  error  signal  represents  the  difference 
between  the  plant  measurement  vector  and  an  identical 
"measurement"  of  the  observer  state  vector.   Specifically 
we  assume  that  the  model  and  the  feedback  error  signal  are 
represented  by  the  observer  state  equations 

i(t)  =  AzCt)  +  G(^(t)  -  Cz(t))         (7D) 

The  feedback  is  provided  by  the  nxm  observer  gain  matrix  G. 
Rewriting  (7D)  in  terms  of  x(t)  and  z^(t)  we  sec  that 
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_z(t)  =  Az(t)  +  GC(xCt)  -  ^Ct)).  (8D) 
Defining  the  error  state  £(t)  as 

eCt)  ^  z(t)  -  xCt)  (9D) 

and  substituting  into  (8D)  yields 

iCt)  =  Az (t)  -  GCe(t) .  (lOD) 

It  is  clear,  therefore,  that  by  using  the  proper  choice  for 
G,  the  error,  e_(t)  ,  can  be  used  to  drive  z_(t)  toward  xC^)  . 

More  precisely,  differentiating  (9D)  with  respect  to 
time  yields 

e(t)  =  i(t)  -  x(t). 

Substituting  (ID)  and  (lOD)  for  x(t)  and  z_(t)  ,  respectively, 
yields 

e(t)  =  AzCt)  -  GCeCt)  -  Ax(t) 

eCt)  =  A(z_(t)  -  x(t))  -  GCe(t) 

e(t)  =  (A-GC)e(t).  (IID) 

That  is 

e(t)  =  Fe_(t)  (12D) 

where  we  define  F  as 

F  =  A  -  GC.  (13D) 

The  solution  to  C12D)  is  given  by 

e(t)  =  e-^e(O)  .  (14D) 

We  see,  therefore,  tliat  if  tlie  eigenvalue?  of  P  all  have 
negative  real  parts,  then  for  a  finite  e(0). 
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lim  e(t)  =  0. 

It  is  shown  in  Re£.  4  that  the  eigenvalues  of  the  ob- 
server can  be  positioned  arbitrarily  if  the  pair  (C,A)  is 
completely  observable.   That  is  if  the  matrix 

has  rank  n.   From  (13D)  we  see  that  the  eigenvalues  of  F  are 
controlled  by  the  selection  of  the  observer  gain  matrix  G. 

The  above  argument  is  easily  extended  to  a  plant  which  is 
driven  by  the  input  u(t) .   The  plant  state  and  output  equa- 
tions become 

x(t)  =  Ax(t)  +  Bu(t)  (15D) 

^(t)  =  Cx(t)  (16D) 

and  the  observer  state  equations  become  ; 

z_Ct)  =  Az(t)  +  G(y(t)  -  Cz(t))  +  Bu(t)(17D) 

where  ju(t)  is  an  £xl  control  vector  and  ^  is  a  known  nxSi 
matrix, 

C.   REDUCED-ORDER  OBSERVERS 

The  full-order  observer  possesses  a  certain  amount  of 
redundancy  since  we  estimate  some  states  which  are  essentially 
being  measured  directly.   We  therefore  desire  to  construct 
an  observer  wliich  estimates  only  those  (n-m)  states  which 
are  not  being  measured  directly.   Wc  start  with  the  state 
and  output  equations 

x(t)  =  Ax(t)  +  Bu(t)  (19D) 
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Y_it)    =  Cx(t) 


C20D) 


We  define  the  new  state  vector  xCt)  so  that  its  first  m 
components  are  the  measurement  vector  ^(t)  and  the  remain- 
ing (n-m)  components  are  the  vector  w(t) .   That  is 


x(t)  ^ 


w(t) 


(21D) 


If  wCt)  is  selected  so  it  is  some  linear  transformation  of 
the  state  x(t),  namely 


w(t)  =  Dx(t) 


(22D) 


Then,  in  view  of  (20D)  and  (22D) ,  (21D)  can  be  written  as 


x(t)  = 


C 
D 


xCt) 


(23D) 


Furthermore,  if  D  is  chosen  such  that 


C 
D 


-1 


exists,  where 


M 


C 
D 


then,  the  state  x(t)  can  be  found  from  x(t) .   Namely 

x(t)  =  M'-^x(t)  . 

Assuming  that  M   exists,  the  task  becomes  one  of  obtain- 
ing an  estimate  for  w(t) .   Writing  the  state  equations  for 
x(t)  we  have 
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That  is 


y_it)    =  CxCt)  =  CAx(t)  +  CBu(t) 
w(t)  =  Dx(t)  =  DAx(t)  +  DBu(t)  . 


x(t)  =  MAx(t)  +  MBu(t) 


(26D) 


Substituting  C25D)  for  x(t)  yields 


Defining 


XCt)  =  MAM'-'-x(t)  +  MBu(t) 


A'  =  MAM'-'-  ,   B'  =  MB 


(27D) 


C28D) 


then  X  =  A'x  +  B'u(t)  (29D) 

which  can  be  written  in  partitioned  matrix  form  as 


w(t) 


A  ' 
t11 

A  ' 
^21 


A  • 
.12 

A  ' 
^22 


Z(t) 
w(t) 


^2 


u(t) 
(30D) 


Therefore 


il(t)  -  A^i^(t)  +  A'2w(t)  +  B^u(t)      (31D) 


and 


w(t)  =  A^;L)^(t)  +  A^2-^^^  "^  B^u(t).     (32D) 
Since  )^(t)  is  a  measured  quantity,  its  time  derivative, 

* 

^(t) ,  can  be  determined.   The  input  u(t)  can  also  be  con- 
sidered as  a  known  quantity.   Therefore,  solving  (31D)  for 
Ai^w(t)  yields 

A'2w(t)  =  ^(t)     -    A[^Z(-^)    -    B'u(t) 
where  the  elements  on  the  right  represent  "known"  quantities 
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The  quantity  Al ^w(t)  can  therefore  be  used  as  the  measurement 
vector  for  a  full-order  observer  for  the  plant  whose  state 
equations  are  given  by  (32D) .   We  therefore  write  the  ob- 
server state  equations  as 

w(t)  =  ^2^^^^    "■  ^21^*^^^  ""  ^2^^^^ 

(33D) 

+  G(y(t)-A^^^(t)-B^uCt)-A^2i(t)) 

where  G  is  an  (n-m)xm  matrix,  w(t)  is  the  estimate  of  w(t) , 

and  where  the  first  three  terms  represent  the  plant  model 

and  the  last  term  represents  the  feedback  term.   Rewriting 
(33D)  yields 

w(t)  =  (^^2'^12^^^^^    ^    CA^rGAi-L)l(t) 

(34D) 
+  (B^-GBpu(t)  +  Gx(t). 

It  has  been  shown  [4]  that  if  (C_,A)  is  observable  then 
(AJ^jA'^)  is  also  observable.   Therefore  G,  the  observer 
gain  matrix,  can  be  selected  so  that  (Aip-GA^2)  ^^^  arbitrary 
eigenvalues . 

Since  differentiation  is  usually  troublesome  to  perform, 
it  is  desirable  to  eliminate  the  term  ^(t)  from  (340) .   To 
accomplish  this  we  subtract  the  term  Gy(t)  from  both  sides 
of  (34D)  to  get 

w(t)-G^(t)  =  (A^2-Mi2^i<;t)  ^  (Ah'^ll^^^^^ 
+  (B^-GBpu(t).  (35D) 

Adding  and  subtracting  (A '  -GA' -) G^(t)  from  the  right  side 
of  (35D)  yields 
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w(t)-Gi(t)  =  (A^2-Mi2Hw(t)-Gx(t)) 


+  (BJ-GBpu(t).  (36D) 


We  now  define 


zCt)  ^  w(t)  -  G^(t) 


Equation  (36D)  can  therefore  be  rewritten  as 

z(t)  =  (A^2-^i2^^'^^^  ^  CA^l-GAi^+(A^2"^i2^^^^<^^^ 

+  (B^-GB|)u(t)          (37D) 
where  w(t)  can  be  computed  from 

w(t)  =  z^(t)  +  Gx(t).  (38D) 
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APPENDIX  E 

SAMPLE  COMPUTER  PROGRAM  FOR  OBSERVER/CONTROLLER  DESIGN 

A  computer  program  is  given  on  the  following  pages  which 
can  be  used  to  design  an  observer/controller  in  the  manner 
prescribed  by  the  algorithm  o£  Section  IID.   It  is  assumed 
that  the  plant  state  and  output  equations  are  given  by 

xCt)  =  AxCt)  +  BuCt) 

Y_(it)    =  Cx(t)  (IE) 

and  that  the  observer  state  and  estimate  equations  are  given 
by 

z_(t)  =  (A-GC)z_(t)  +  G>^Ct)  +  Bu(t).  (2E) 
The  control  law  is  given  by 

u(t)  =  Kz(t)  (3E) 

and  the  performance  measure  to  be  minimized  is  given  by 


J  =  /  t^[x'^(t)gx(t)  +  e^Ct)Re(t)]dt.  (4E) 


The  following  symbols  are  explained  in  the  order  of  their 
appearance  in  the  input  section  of  the  program: 

BETA  -  The  initial  step  length  for  the  minimizing  search 
routine  DDIRCT. 

EPS  -  The  minimum  step  length  allowed. 

D  -  The  minimum  value  that  the  sum  of  tlie  eigcn,values 
(trace)  of  the  closed-loop  system  is  allowed  to  take  on. 
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ALPH  -  The  weighting  multiplier  which  is  used  to  assign 
a  cost  to  any  violation  of  the  trace  constraint. 

N  -  The  order  of  the  matrix  A  in  (IE). 

JY  -  The  maximum  number  of  searches  the  minimizing  search 
routine  is  allowed  to  perform. 

NR  -  The  time  multiplier  exponent  r  in  the  performance 
measure  (4E) . 

M  -  The  row  dimension  of  C  in  (IE) .   For  full-order 
observer/controller  design  M  is  set  to  zero. 

A  -  The  matrix  A  in  (IE)  . 

Q  -  The  matrix  qi  in  (4E)  . 

BB  -  The  matrix  B  in  (IE). 

GG  -  A  vector  whose  elements  are  given  by  the  initial 
estimate  of  G_  and  K  in  the  form 

G  " 
K 

SS  -  The  R  matrix  in  (4E) 
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